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Abstract. A global analysis is presented of solutions for Laplace's equation on three- 
dimensional Euclidean space in one of the most general orthogonal asymmetric confocal 
cyclidic coordinate systems which admit solutions through separation of variables. We 
refer to this coordinate system as five-cyclide coordinates since the coordinate surfaces 
are given by two cyclides of genus zero which represent the inversion at the unit sphere 
of each other, a cyclide of genus one, and two disconnected cyclides of genus zero. This 
coordinate system is obtained by stereographic projection of sphero-conal coordinates 
on four-dimensional Euclidean space. The harmonics in this coordinate system arc 
given by products of solutions of second-order Fuchsian ordinary differential equations 
with five elementary singularities. The Dirichlet problem for the global harmonics in 
this coordinate system is solved using multiparameter spectral theory in the regions 
bounded by the asymmetric confocal cyclidic coordinate surfaces. 
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1. Introduction 

In 1894 Maxime Bocher's book "Ueber die Reihenentwicklungen der Potentialtheorie" 
was published [2]. It took its origin from lectures given by Felix Klein in Gottingen 
(see for instance IE])- In Bocher's book, the author gives a list of 17 inequivalent 
coordinate systems in three dimensions in which the Laplace equation admits separated 
solutions of the form 



where the modulation factor R(x,y,z) [HI p. 519] is a known and fixed function, and 
Si, s 2 , s 3 are curvilinear coordinates of x, y, z. The functions W\, w 2 , w 3 are solutions of 
second order ordinary differential equations. The symmetry group of Laplace's equation 
is the conformal group and equivalence between various separable coordinate systems 
is established by the existence of a conformal transformation which maps one separable 
coordinate system to another. 

In general, the coordinate surfaces (called confocal cyclides) are given by the zero 
sets of polynomials in x, y, z of degree at most 4 which can be broken up into several 
different subclasses. For instance, eleven of these coordinate systems have coordinate 
surfaces which are given by confocal quadrics (9j Systems 1-11 on p. 164], nine are 
rotationally-invariant [9, Systems 2,5-8 on p. 164 and 14-17 on p. 210], four are 
cylindrical [9], Systems 1-4] , and five which are the most general, are of the asymmetric 
type namely, confocal ellipsoidal, paraboloidal, sphero-conal, and two cyclidic coordinate 
systems |HJ Systems 9,10,11 on p. 164 and Systems 12,13 on p. 210]). Bocher [2] showed 
how to solve the Dirichlet problem for harmonic functions on regions bounded by such 
confocal cyclides. However, it is stated repeatedly in Bocher's book that the presentation 
lacked convergence proofs, for instance, this is mentioned in the preface written by Felix 
Klein. 

It is the purpose of this paper to supply the missing proofs for one of the asymmetric 
cyclidic coordinate systems which is listed as number 12 in Miller's list P, page 210] (see 
also Boyer, Kalnins & Miller (1976) [3, Table 2] and Boyer, Kalnins & Miller (1978) [I] for 
a more general setting). For lack of a better name we call it 5-cyclide coordinates. This 
asymmetric orthogonal curvilinear coordinate system has coordinates s, G R (i = 1, 2, 3) 
with Si in (a ,cii), (ai,a 2 ) or (a 2 ,a 3 ), respectively where a < ai < a 2 < a 3 are 
given numbers. This coordinate system is described by coordinate surfaces Sj = const 
which are five compact cyclides. The surfaces si = const for sx G (a ,ai) are two 
cyclides of genus zero representing inversions at the unit sphere of each other. The 
surface s 2 = const for s 2 G (ai,a 2 ) represents a ring cyclide of genus one and the 
surfaces s 3 = const for s 3 G (a 2 , a 3 ) represent two disconnected cyclides of genus 
zero with reflection symmetry about the x, y-plane. The asymptotic behavior of this 
coordinate system as the size of these compact cyclides increases without limit is 6- 
sphere coordinates (see Moon & Spencer (1961) pm p. 122]), the inversion of Cartesian 
coordinates. 




(1) 
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In our notation the coordinate surfaces of this system are given by the variety 

(x 2 + y 2 + z 2 - l) 2 Ax 2 Ay 2 Az 2 , , 

+ + + = 0, (2) 

s — do s — a% s — 02 s — as 

where s = Sj is either in (a , fli), (oi,o 2 ) or (02,03) respectively See Figure fl](a,b) for 
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(a) s 1= 0.7, s 2 =L01, 53=2.01 


(b) ^=0.01, 52=1.5, s 3 =2.5 



Figure 1. Surfaces Si.2,3 = const for — i, where only the component of the cyclide 
si = const inside the ball x 2 + y 2 + z 2 < 1 is shown. 



a graphical illustration of these triply-orthogonal coordinate surfaces, where we have 
selected one of the confocal cyclides for si = const. This is a very general coordinate 
system containing the parameters 0,0,01,02,03 which generates many other coordinate 
systems by limiting processes. Since the book by Bocher is quite old and uses very 
geometrical methods, we will present our results independently of Bocher's book. We 
supply convergence proofs based on general multiparameter spectral theory [H [13] which 
was created with such applications in mind. As far as we know this general theory has 
never before been applied to the Dirichlet problems considered by Bocher. 

We start with the observation that 5-cyclide coordinates are the stereographic 
image of sphero-conal coordinates in four dimensions (or, expressed in another way, of 
ellipsoidal coordinates on the hypersphere S 3 ). We take the sphero-conal coordinate 
system as known but we present the needed facts in Section [2] The well-known 
stereographic projection is dealt with in Section [3] which also explains the appearance 
of the factor R in Q. The 5-cyclide coordinate system is introduced in Section [5J The 
solution of the Dirichlet problem on regions bounded by surfaces ^ with s E (ai, 02) is 
presented in Section [6j The preceding Section [5] provides the needed convergence proofs 
based on multiparameter spectral theory. The remaining sections treat the Dirichlet 
problem on regions bounded by the surfaces §2§ when s G (oi,a 2 ) (ring cyclides) and 

s e (o 2 ,o 3 ). 
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2. Sphero-conal coordinates on R fc+1 

Let k 6 N. In order to introduce sphero-conal coordinates on R fc+1 , fix real numbers 

do < CL\ < d2 < ■ ■ ■ < Ok- (3) 

Let (x , Xi, . . . , Xk) be in the positive cone of R fe+1 

x > 0,...,x k > 0. (4) 
Its sphero-conal coordinates r, s±, . . . , Sf. are determined in the intervals 

r > 0, aj_! < Si < a,;, i — 1, . . . , k (5) 
by the equations 

k 

r 2 = ]>>? (6) 

3=0 

and 

k 2 

~z = fori = l,...,Jfe. (7) 
j=o s * _ 

The latter equation determines si, s 2 , ■ ■ ■ , s k as the zeros of a polynomial of degree k 
with coefficients which are polynomials iniQ,...,!^. 

In this way we obtain a bijective (real-)analytic map from the positive cone in R fe+1 
to the set of points (r, si, . . . , satisfying The inverse map is found by solving a 
linear system. It is also analytic, and it is given by 

2 _ 2 rL=i( g » ~ Q j) /o\ 

Sphero-conal coordinates are orthogonal, and its scale factors (metric coefficients) are 
given by H r = 1, and 

Hi = \ V = -] r2 U T liSl - , i = l,2,...,Jfe. (9) 

1 4 £z (* - 4 n -=o(^ - «i) 



Consider the Laplace equation 

dx 



^ = Eg = o do) 

for a function C/(s ,^i, . . . , x^). Using ([£]) we transform this equation to sphero-conal 
coordinates, and then we apply the method of separation of variables [12] 



U (x , xi, . . . , x k ) = w (r)w 1 (s 1 )w 2 (s 2 ) . . . w k (s k ). (11) 
For the variable r we obtain the Euler equation 
// A.' i 4Ao 

w o + ~ w o H 5" w o = (12) 
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while for each of the variables s\, s 2 , 



j=0 

More precisely, if A , 



w 



1 k 

IE 



s — a* 
o 



Sk we obtain the Fuchsian equation 

fc-i 

Ev- 1 - 



1 



-w 



+ 



0. 



(13) 



, Afc_i are any given numbers (separation constants), and if 



u>o(r), r > 0, solves (12) and Wi 



-i < Si < a,i, solve (1131 for each z = 1, 



k, then 



£/ defined by (11) solves (10) in the positive cone of R fc+1 rt4b. 



Equation ( 13 ) has only regular points except for k + 2 regular singular points at 

oo. The exponents at each finite singularity s 



, cifc and s 



dj are and 



s — cto, &i , 

|. Therefore, for each choice of parameters Aq, ■ ■ ■ , Afe_i, there is a nontrivial analytic 



solution at s = aj and another one of the form w(s) 
at aj. If v, /i denote the exponents at s = oo then 



a 



a/a 



An 



The polynomial Xlto 1 AjS fc_1_ * appearing in (|l3j) is known as van Vleck polynomial. If 
k = 1 then (13) is the hypergeometric differential equation (up to a linear substitution). 
If k = 2 then (13) is the Heun equation. We will use this equation for k = 3. According 
to Miller (1977) p. 209] (see also [3, p. 71]) in reference to the k = 3 case, "Very little 
is known about the solutions." 



yU + v 



k 



v(s), where v is analytic 



(14) 



3. Stereographic projection 

We consider the stereographic projection P : § 3 \ {(1, 0, 0, 0)} — > R 3 given by 

P(x ,Xi,x 2 ,x 3 ) = - — — {x 1 ,x 2 ,x 3 ). 
l-x 

The inverse map is 

P- I (x, y, z) = — 2 \ (x 2 + y 2 + z 2 -l, 2x, 2y, 2z). 

x z + y z + z z + 1 

We extend P^ 1 to a bijective map 

Q : (0, oo) x R 3 -> R 4 \ {(x , 0, 0, 0) : x > 0} 

by defining 

Q(r,x,y,z) := rP _1 (x, y, z). 

If we set (xo, Xi, X2, x 3 ) = Q(r, x, y, z) we may consider r, x, y, z as curvilinear coordinates 
on R 4 with Cartesian coordinates x ,Xi,x 2 ,x 3 . We note that x^ + x\ + x 2 + x\ = r 2 so 
r is just the distance between (xo, Xi, x 2 , Xg) and the origin. Moreover, (x,y,z) is the 
stereographic projection of the point (xo/r,Xi/r,x 2 /r,x 3 /r) G § 3 . It is easy to check 
that the coordinate system is orthogonal and scale factors are 

h r = 1, h x = h y = h z = 2rh, where h := —„ 5— — r -. 

x 2 + y 2 + z 2 + 1 
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Let U(xq, Xi,X2,X3) = V(r,x,y, z). Then 

AU = {(2rhV x ) x + (2rhV y ) y + (2rhV g ) z + (8r 3 h 3 V r ) r ) . (15) 

Suppose that U is homogeneous of degree a: 

U (txo, txi,tx 2 ,tx 3 ) = t a U(x , Xi, x 2 , x 3 ), t > 0. 
Then V can be written in the form 

V(r,x,y,z) = r a w(x,y,z), 



and (15) implies 

r a-2 

AU = W (( hw *)* + ( hw y^y + ( hw ^ z + 4a ^ a + 2 ) h3w ) ■ ( 16 ) 

We now introduce the function 

u(x, y, z) = w(x, y, z)(x 2 + y 2 + z 2 + 1) -1 ^ 2 . 



Then a direct calculation changes (16) to 

AU = (u xx + u yy + u zz + (3 + Aa{a + 2))h 2 u) . (17) 

If 3 + 4a (a + 2) = then U is harmonic if and only if u is harmonic. Noting that 
3 + 4a (a + 2) = (2a + l)(2a + 3), we obtain the following theorem. 

Theorem 3.1. Let D be an open subset of S> 3 not containing (1,0,0,0), let E = 
{(rxo, rxx,rx2, rx%) : r > 0, (xq, Xi, x 2 , £3) G D}, and let F = P(D) be the stereographic 
image of D. Let the function U : E — >■ R be homogeneous of degree — | or — |, and 
let w : F R satisfy U = w o P on D. Then U is harmonic on E if and only if 
w(x, y, z){x 2 + y 2 + z 2 + l)^ 1 / 2 is harmonic on F. 

4. Five-cyclide coordinate system on R 3 

We introduce sphero-conal coordinates 

r > 0, ao < si < ai < s 2 < «2 < S3 < 03, 

on R 4 as explained in Section [2] with k = 3. Then si,s 2 ,ss form a coordinate system 
for the intersection of the hyper sphere S> 3 with the positive cone in R 4 . Using the 
stereographic projection P from Section [3] we project these coordinates to R 3 . We 
obtain a coordinate system for the set 

T — {(x, y, z):x,y,z> 0, x 2 + y 2 + z 2 > 1}. (18) 



Explicitly, 



where 



X\ x 2 2 3 

x = - , y=- , z = - , (19) 

1 — X 1—^0 1 — 



2= UWi %) J = ,l,2,3. (20) 

n&i^w ~ a j) 
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Conversely, the coordinates s\, S2, S3 of a point (x, y,z) G T are the solutions of 

0. (21) 



(x 2 + y 2 + z 2 - l) 2 ^ Ax 2 



Ay 2 Az 2 



s — ao s — ai s — CL2 s — a 3 

Since sphero-conal coordinates are orthogonal and the stereographic projection preserves 
angles, 5-cyclide coordinates are orthogonal, too. This is the twelfth coordinate system 
in Miller (1977) j9j page 210]. Miller uses a slightly different notation: a — 0, a± — 1, 
02 = b, a 3 = a, and si = p, S2 = v, S3 = \i. Also, x, z are interchanged. 

In order to calculate the scale factors for the 5-cyclide coordinate system we proceed 
as follows. We start with 



dx 
dsi 



dxi 



+ 



dXr 



1 - x dsi ' (1 - x ) 2 dsi ' 
and similar formulas for the derivatives of y and z. Then using 



Xq ~\~ X^ ~\~ X2 ~\~ x^ 



dxo dx\ 8x2 8x3 

1" Xi— h X 2 ^. h X3- 



dsi 



dsi 



a short calculation gives 



dx dx dy dy <Jz <)z 



dsj dsj dsi dsj 



dsi dsj 



dsi 



dxo dxf 



dsi 



Xq) 2 ' dsi ds^ 



1=0 



This confirms that 5-cyclide coordinates are orthogonal and from ^ we obtain the 
squares of their scale factors 



hi 



1 

16 



(P 2 " 1] 



Si 



a 



+ 



Ax 1 



+ 



% 2 



+ 



Az 2 



Si - ai) 2 (si - a 2 ) 2 (si - a 3 



where p 2 = x 2 + y 2 + z 2 , or, equivalently 
1 

16 
1 

16 
1 

16 



hi 



h 2 2 



hi 



:(P 2 + 1) 2 
;(P 2 + 1) 2 

V + i) 2 



(«8 


- s x ){s 2 


-si) 


(s 1 - a )(a 1 


- si)(a 2 


- si)(a 3 - s x ) 




- si)(s 3 


-S2) 


(s 2 - a )(s 2 


- ai)(a 2 


- s 2 )(a 3 - s 2 ) 


(S3 


- si)(s 3 


-82) 


(S3 - ao)(s3 


- ai)(s 3 


~ a 2 )(a 3 - S3) 



(22) 

(23) 
(24) 
(25) 



We find harmonic functions by separation of variables in 5-cyclide coordinates as 
follows. 



Theorem 4.1. Let w\ : (ao,ai) 
of the Fuchsian equation 
3 



C, W2 '■ (01,0-2) — > C, 103 : (a 2 ,a 3 ) — > C be solutions 



I* 

3=0 



s — a,- 



w" + l~y~] — - — w' + ( ^-s 2 + Xis + A 2 ) w 
2 ^ s — a,- \ 16 / 



0. 



where \±,\2 are given (separation) constants. Then the function 
u(x, y, z) = (x 2 + y 2 + z 2 + l)" 1/2 wi(s 1)^2(^2)^3(^3) 



(26) 



(27) 



is a harmonic function on the set (18). 
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Proof. Using sphero-conal coordinates r, s\, S2, S3 on R 4 , we define a function U in the 
positive cone of R 4 by 

U(x ,x 1 ,x 2 ,x 3 ) = r~ 1/2 w 1 (s 1 )w 2 (s 2 )w 3 (s 3 ). 

The function r -1 / 2 is a solution of (12) when k = 3, Xq = j^. The results from Section^ 
imply that U is harmonic, and, of course, U is homogeneous of degree — |. The function 
u> defined on the set ( 18 ) by U = w o P is given in 5-cyclide coordinates by 



w(x,y,z) = w 1 (s 1 )w 2 (s 2 )w 3 (s 3 ). 



Therefore, Theorem 3.1 gives the statement of the theorem. 



□ 



Equation (26) has five regular singularities at s = ao,ai,a 2 ,a 3 ,oo. The exponents 
at the finite singularities are and |. Using (14), we find that the exponents at infinity 
are j and |. So all five singularities are elementary in the sense of Ince [5]. Equation 
(26) is one of the standard equations in the classification of Ince [5J page 500]. 

We define the 5-cyclide coordinates si,S2, s 3 for an arbitrary point (x, y,z) G R 3 
as the zeros si < s 2 < s 3 of the cubic equation 



3 

Y[(s - aj 

j=0 

For example, Sj (0, 0, 0) 



(x 2 + y 2 + z 2 - l) 2 4x 2 
s — do s — at 



+ 



% 2 



a 2 



+ 



4z 2 



a 3 



0. (2£ 



Qj for j = 1,2,3. Each function Sj : R 3 — > [flj_i,aj 



is 

continuous. We observe that, in general, there are 16 different points in R 3 which have 
the same coordinates s%,s 2 ,s 3 . If (x, y,z) is one of these points the other ones are 
obtained by applying the group generated by inversion at S 2 

a (x,y,z) = p~ 2 (x,y,z) (29) 

and reflections at the coordinate planes 

<ri(x,y,z) = (-x,y,z), <r 2 (x,y, z) = (x,-y,z), a 3 (x,y,z) = (x,y,-z).(30) 



It is of interest to determine the sets where Sj = a,_i or Sj 



si 

S3 
S3 



a iff x 2 + y 2 + z 2 
and 



a\ iff x 

a\ iff x 

a 2 iff y 

a 2 iff y 
a 3 iff z 



and 

and 

and 
0. 



(P 2 


-I) 2 


+ 


Ay 2 


+ 


Az 2 


a 1 


- ao 




d\ — a 2 


d\ — a 3 


(P 2 


-I) 2 


+ 




+ 


Az 2 


a 1 


- ao 


a\ — a 2 




ai — a 3 


(P 2 


-I) 2 


+ 


Ax 2 


+ 


Az 2 


a 2 


- ao 


a 2 — ai 


d 2 — a 3 


(P 2 


-I) 2 


+ 


Ax 2 


+ 


Az 2 


a 2 


- ao 


a 2 — a\ 


a 2 — a 3 



dj. We obtain 



>0, 
<0, 

>o, 
< 0, 



We define the sets (consisting each of two closed curves) 



A x := {(x, y, z) G R 3 '■ S\ = s 2 = di} 



(31) 
(32) 

(33) 

(34) 

(35) 
(36) 

(37) 



{{x,y,z) : x = 0, 



(P 2 



CL\ — Oq 



+ 



Ay 2 



Az 2 



+ 

d\ — CL 2 CLi — 03 



0}, 
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see Figure |2j and 

M ■= {{x, y, z) E R 3 : s 2 = s 3 = a 2 } 



(38) 



{{x,y,z) :y = 0, 



(P 2 " If 



Ax 2 Az 2 
+ 



0}, 



0,2 — Oq 0-2 — 0\ 2 — G<3 

see Figure [3] Clearly, sj is analytic at all points (x, y, z) at which sj is a simple 




Figure 2. The set Ai in the y,z-plane for ai = i. 



z 




Figure 3. The set A2 in the x, z-plane for a; = i. 
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zero of the cubic equation (28). Therefore, si is analytic on R \A\, S2 is analytic on 



R 3 \ (Ai U A 2 ), and s 3 is analytic on R 3 \ A 2 . 



We may use (27) to define u(x,y,z) for all (x,y } z) G R . Since the solutions 
Wi,W2,W3 of (26) have limits at the end points of their intervals of definition (because 
the exponents are and | there), we see that u is a continuous functions on R 3 . The 
function [x 2 + y 2 + z 2 + l) 1 / 2 w(a;, y, z) is invariant under Cj, % — 0, 1,2,3. In general, 
u is harmonic only away from the coordinate planes and the unit sphere. In fact, we 
observe that u is a bounded function which converges to at infinity, so, by Liouville's 
theorem, u cannot be harmonic on all of R 3 unless it is identically zero. 

5. First two-parameter Sturm-Liouville problem 



We consider equation (26) on the intervals (01,02) and (02,03) and write it in formally 



self-adjoint form. Setting 

1 /2 

u(s) := \(s - a )(s - ai)(s - o 2 )(s - a 3 )| , (39) 
we obtain two Sturm-Liouville equations involving two parameters 

(ou(s 2 )w' 2 )' + — j—r ( ^sl + Xis 2 +A 2 J w 2 = 0, a 1 <s 2 <a 2 , (40) 

(u)(s 3 )w' 3 )' ^— ( ^-s 2 3 + A1S3 + A 2 J w 3 = 0, a 2 < s 3 < a 3 . (41) 

u(s 3 ) V 16 / 



In (40) W2 is a function of S2 and the derivatives are taken with respect to s 2 . In (41) 
W3 is a function of S3 and the derivatives are taken with respect to S3. We simplify the 
equations by substituting tj = Q(sj), Uj(tj) = Wj(sj), where f2(s) is the elliptic integral 
(see for instance [S]) 

m--=[ S ^ry (42) 

This is an increasing absolutely continuous function Q : [00,03] — > [0,63], where 



bj := 0(aj). Let <p : [0,6 3 ] — > [oo,o 3 ] be the inverse function of Q. Then (40), (41) 
become 

«2+ ^W*2)} 2 + V(i 2 )U 2 j« 2 = 0, h<t 2 <b 2 , (43) 

u'i - (^{4>(h)} 2 + Ai0(t 3 ) + A 2 J % = 0, 6 2 < * 3 < 63. (44) 
We add the boundary conditions 

4(60 = u' 2 (b2) = 1/3(62) = 4(6 3 ) = 0. (45) 



Differential equations (43), M4J) together with boundary conditions (45) pose a two- 
parameter Sturm-Liouville eigenvalue problem. For the theory of such multiparameter 
problems we refer to the books [HUH] and the references therein. A pair (Ai, A2) is called 
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an eigenvalue if there exist (nontrivial) eigenf unctions u^it-z) and ^3(^3) which satisfy 
(43), (44), (45). The two-parameter problem is right- definite in the sense that 

1 



0(* 2 ) 
-<f>{h) 



■1 



0(t 3 ) - 0(t 2 ) > for 61 < t 2 < b 2 < t 3 < 63. 



However, this determinant is not positive on the closed rectangle [61,62] x [62,63]. This 
lack of uniform right- definiteness make some proofs in this section a little longer than 
they would be otherwise. 

We have the following Klein oscillation theorem; see [TJ Theorem 5.5.1]. 

Theorem 5.1. For every n = (n 2 ,n 3 ) G Nq ; there exists a uniquely determined 
eigenvalue (Ai jn , A2, n ) £ R 2 admitting an eigenfunction u 2 with exactly n 2 zeros in 
(61,62) and an eigenfunction u% with exactly n% zeros in (62,63). 

We state a result on the distribution of eigenvalues; compare with [TJ Chapter 8]. 

Theorem 5.2. There are positive constants A\, A 2 , A 3 , A4 such that, for all n G Nq ; 

- A x {n\ + n\ + 1) < Ai, n < -A 2 {n\ + n\) + A 3 , (46) 

|A 2j „| < Ai(n 2 2 + n 2 3 + l). (47) 



Proof. If a differential equation u" + q(t)u = with continuous q : [a, b] — y R admits a 
solution u satisfying u'(a) = u'(b) = and having exactly m zeros in (a, 6), then there 

2 2 

is t e (a, 6) such that g(t) = 5_^. 2 . This is shown by comparing with the eigenvalue 



problem u" + \u = 0, u'(a) = u'(b) = 0. Applying this fact, we find t 2 G (61,62) and 
t 3 G (6 2 , 63) such that 



3 

16 

3 

16 



{(p(t 2 )} 2 + XMh) + \ 2 

{0(t3)} 2 + A!0(t3) + A2 



2^2 



7T n 



(62 -61) 2 ' 



(63 - 6 2 ) 



(48) 
(49) 



where we abbreviated Xj = Xj n . By subtracting (48) from (49), we obtain 



16 



(W3)} 2 - W2)} 2 ) + Ai(0(t 3 ) - 0(t 2 )) 



2^2 



7r n 



2^2 



7r n 



(62 - 61 



(63 - 6 2 ) 



< 0. 



Dividing by ^(£3) — 4>(t 2 ) and using < ^(£3) — 0(^2) < 03 — ai, we obtain the second 
inequality in (46). 



To prove the first inequality in (46), suppose that Ai < — |a 3 . Then the van Vleck 
polynomial 

3 



Q{s) 



16 



s 2 + Ais + A 2 



(50) 



satisfies Q'(s) = |s + Ai < for s < a 3 . Let c G (61,62) be determined by 
0(c) = |(ai + a 2 ). If Q(fl2) > then, for t G [61, c], 

3 



Q(0(t))>Q(^K + a 2 ))>^(a 2 



-Ai 



:«3 
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By Sturm's comparison theorem applied to equation (43), we get 
(c - b x f{a 2 - a x ) (-Xx - -a 3 ) < 4vr 2 (n 2 + l) 2 , 



which gives the desired inequality. If Q{a 2 ) < we argue similarly working with (44) 
instead. 

Finally, g7)) follows from Q and (g. □ 

Let M2, n and u 3yXl denote real-valued eigenfunctions corresponding to the eigenvalue 
(^i,m -^2,n)- It is known [U section 3.5] (and easy to prove) that the system of products 
u 2 , n (t 2 )u 3}Il (t 3 ) , n G Nq, is orthogonal in the Hilbert space Hi consisting of measurable 
functions / : (61,62) x (62,63) — > C satisfying 

(<f)(t 3 ) - <f>{t 2 )) \f{t 2 , t 3 )\ 2 dt 2 dt 3 < 00 



62 

with inner product 

fb 3 rb 2 



b 2 



((j)(t 3 ) - (j)(t 2 ))f(t 2 , t 3 )g(t 2 ,t 3 ) dt 2 dt 3 . 



We normalize the eigenfunctions so that 

-62 



b 2 Jb\ 



(<f>(t 3 ) ~ <f>{t 2 )) {^2,n(t 2 )} 2 {U 3 , n (t 3 )f dt 2 dt 3 = 1. 



(51) 



We have the following completeness theorem; see [13, Theorem 6.8.3]. 
Theorem 5.3. The double sequence of functions 

U2, a {h)u 3>a {t 3 ), n G Nq, 
forms an orthonormal basis in the Hilbert space H\. 



The normalization (51) leads to a bound on the values of eigenfunctions. 

Theorem 5.4. There is a constant B > such that, for all n G Nq and all t 2 G [61, 62], 
h e [62,63], 

\u 2 , n (t 2 )u 3}ll (t 3 )\ < B{n\ + n\ + 1). 



Proof. We abbreviate Uj = Uj jX1 , Xj = Xj jU . Condition (51) is a normalization for 
the product ^2(^2)^3(^3) but not for each factor separately, so we may assume that, 
additionally, 

(52) 



[ 2 {u 2 (t 2 )} 2 dt = 1. 

Jbi 



Now ([511 , (|52j) imply that 

rb 3 



(0(t3)-0(6 2 )){n 3 (t3)r dt 3 <l. 



(53) 



b-2 
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We multiply equations (43), (44) by u 2 and w 3 , respectively, and integrate by parts 
to obtain 



fb 2 






12 




u 2 


Ai 










12 






lb. 





3 

16 



l>2 



<?u\ + Ai 



b-2 



by 



a 



2- 



3 

16 



6 3 



- Ai 



l>2 



63 



*>u 3 - A 2 



b-2 



b?, 



(54) 



(55) 



6 2 



It follows from (52), (54) and Theorem 5.2 that there is a constant B\ > such that, 
for all n G Nq, 



pb 2 

/ U? < #1(7^ + ^3 + 1). 



(56) 



Unfortunately, we cannot argue the same way for U3 because we do not have an upper 

fb'i 

equations. Then, noting (BTj), we find 

fb 2 r b '-i 



bound for J b b3 u\. Instead, we multiply (54) by J u\ and (55) by J u\ and add the 



12 
u 2 



62 



/•b 2 rbz g 

+ / u \\ % 2 < -Ai + - max 10(f) I . 



Using Theorem 5.2 and (52), we find a constant B 2 > such that, for all n G Nq, 



b 3 



u 3 2 < B 2 (nj + n% + l). 



(57) 



b 2 



We apply the following Lemma 5J)] (noting (52), (53), (56), (57)) and obtain the desired 
result. □ 

Lemma 5.5. Let u : [a, b] — > R be a continuously differentiate function, and let 
a < c < d < b. Then, for all t G [a, b], 

fd rb 

(d-c) \u{t)\ 2 < 2 / \u(r)\ 2 dr + 2(b - a)(d - c) / \u\r)\ 2 dr. 



Proof. For s,t G [a, b] we have 



\u(t) — u(s)\ 



u{r) dr 



< \t-s\ 1/2 I / \u'(r)\ 2 dr 



1/2 



\u(t)\ 2 < 2\u(s)\ 2 + 2\t - s\ / \u'{r)\ 2 dr. 



This implies 



We integrate from s = c to s = d and obtain the desired inequality. 
Let Mi n be the solution of 



16 



W1)} 2 + XtMt) + A 2 , n )u 1 = 0, 6 < tx < h, 



□ 



(58) 



determined by the initial conditions 

u 1 (b 1 ) = l, u'M = 0. 
The following estimate on w ln will be useful in Section [6j 
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Theorem 5.6. We have «i >n (£i) > for all t\ G [&o>&i]- If = b < c± < C2 < b\, then 
there are constants C > and < r < 1 such that, for all n G Nq and ti G [c2, &i], 

^l,n(^l) < ^ r n 2 +ri3 
Ml,n(ci) ~~ 

Proof. We abbreviate Mi = wi in and Aj = Aj in . By definition, u\ satisfies the differential 
equation 

u'[ = Q((f>(t 1 ))ui, h G [&oAL 



where Q is given by (50). According to (48), (49), there are S2 G (cti, 02) and S3 G (a 2 , 03) 
such that 



n 2 nn „ , . Ti 2 n\ 



(62 - &1) 2 ' " VJ/ (63-6 



2 J 



\2' 



If s < S2 then Q(s) > L(s), where L(s) is the linear function with L(sj) = Q(sj), 
j = 2, 3. It follows that Q(s) > for s G [a , a\\ and 

> C x (nz + n 3 ) 2 for * x G [6 , c 2 ], (59) 

where C\ is a positive constant independent of n. We now apply the following Lemma 



5.7 to complete the proof. Note that the interval [a, b] in the lemma is [ci,&i] but with 



the end points interchanged. □ 

Lemma 5.7. Let u : [a, b] — > R be a solution of the differential equation 

u"{t) = q(t)u(t), te[a,b], 

determined by the initial conditions u(a) = 1, u'(a) = 0, where q : [a, b] — > R is a 
continuous function. Suppose that q(t) > on [a,b] and q(t) > A 2 on [c,b] for some 
A > and c G [a, b). Then u(t) > for all t G [a, b], and 

^2 > 1 a(6-c) j allte[0,c\. 
u{t) - 2 J L ' J 

Proof. Since g(t) > 0, w(t) > and u'(t) > for t G [a, 6]. The function z = -u'/w 
satisfies the Riccati equation 

z' + Z 2 = q(t), 
and the initial condition z(a) — 0. It follows that 

z(t) > Atanh(A(t-c)) for t G [c, b]. 
Integrating from t = c to t = b gives 



In 



u{b) 



> lncoshA(6-c) > In -e A(b " c) 



u(c) ~ ~ 2 

which yields the claim since u is nondecreasing. □ 
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We now introduce a systematic notation for our eigenvalues and eigenfunctions. 
First of all, we note that the results of this section remain valid for other sets of 
boundary conditions. We will need eight sets of boundary conditions labeled by 
P — (PiiP2,P3) G {0,1} 3 - These boundary conditions are 

w / 2 (6 a ) = ifpi = 0, w 2 (6i) = if Pi = 1, 

u' 2 {b 2 ) = u' 3 {b 2 ) = if p 2 = 0, u 2 (b 2 ) = u 3 (b 2 ) = if p 2 = 1, (60) 
u' 3 {b 3 ) = ifp 3 = 0, u 3 (b 3 ) = if Pa = I- 

The initial conditions for u\ are 

ui(6i) = 1,^(60 = ifpi = 0, ui(6 1 ) = 0,«i(6 1 ) = l ifpi = l. (61) 

We denote the corresponding eigenvalues by (A^„ )P , A 2 |„_ p ). For the notation of 
eigenfunctions we return to the Sf variable connected to U by ti = il(si). The 
eigenfunctions will be denoted by E^ p (si) = u i)Tl (ti), i = 1,2,3. The superscript 
(1) is used to distinguish from eigenvalues and eigenfunctions introduced in Sections 7 
and 9. The subscript n = (n 2 , n 3 ) indicates the number of zeros of ^ 2 n P ( s 2), ^3iip( s 3) 
in (ai, a 2 ), (a 2 , as), respectively. The subscript p indicates the boundary conditions used 
to determine eigenvalues and eigenfunctions. By using the letter E for eigenfunctions 
we follow Bocher [2J. In our notation we suppressed the dependence of eigenvalues and 
eigenfunctions on ao, cli, a 2 , a 3 . 



Summarizing, for i — 1, 2, 3, is a solution of (|26|) on (aj_i, a^) with (Ai, A 2 ) 



(Ainp5A 2 n P ). The solution E\^ p (si) has exponent \p x at a± and it has no zeros in 

2,n,pl 

?7 3,n,p^ 



(a , ai). The solution -E 2 1 n iP (s 2 ) has exponent |pi at a 1; exponent |p 2 at a 2 , and its has 
n 2 zeros in (a 1; a 2 ). The solution E^l v (s 3 ) has exponent |p 2 at a 2 , exponent at a 3 , 
and it has n 3 zeros in (a 2 ,a 3 ). 

6. First Dirichlet problem 



Consider the coordinate surface (21) for fixed s = d\ G (a ,ai). See Figure |4j for 

a graphical depiction of the shape of this surface. Let (x r , y', z') G S> 2 . The ray 
(x, y, z) = t(x', y', z'), t > 0, intersects the surface if 

f' 2 - 1 ' 2 -^ (62) 



d\ — Gq 



where 



Ax' 2 Ay' 2 Az' 2 

+ ; + r- > 0. 



a\ — d\ a 2 — d\ a 3 — d\ 
Equation (62 ) has two positive solutions t = ti, t 2 such that < t\ < 1 < t 2 and t\t 2 = 1. 
Therefore, the coordinate surface s± = d\ consists of two disjoint closed surfaces of genus 
0. One lies inside the unit ball -Bi(O) centered at the origin and the other one is the 



image of it under the inversion (29). Let D\ be the region interior to the first surface, 
that is, 

D 1 = {(x,y,z)EB 1 (0):s 1 >d 1 }, (63) 
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(a) s 1= 0.99 




(b) s 1= 0.5 




(c) s 1= 




(d) S!=0.5 




(e) 5i=0.9 




(f) si=0.99 



Figure 4. Coordinate surfaces si = const for a, = z with (a), (b) inside -Bi(O), (d), 
(e), (f) outside £?i(0), and (c) the unit sphere. 



or, equivalently, 



D 1 = {{x,y,z)eB 1 {0): 



(P 2 



r 



+ 



% 2 



+ 



4z 2 



>0}. 



We showed that D\ is star-shaped with respect to the origin. We now solve the Dirichlet 
problem for harmonic functions in D\ by the method of separation of variables. 

Let p = (pi,P2)Ps) ^ {0, l} 3 and n = (n 2 ,n 3 ) G Nq. Using the functions E^ p 
introduced in Section [5] we define the internal 5-cyclidic harmonic of the first kind 

G<$ p {x,y,z) = (x^^^l)- 1 / 2 ^ (64) 

for x,y,z G #l(0) with x,y,z > 0. We extend this function to -Bi(O) as a function of 
parity p. We call a function / of parity p if 



f(<Ti(x,y,z)) = (-l) p *f(x,y,z), 



for i 



1,2,3 



(65) 



using the reflections ( 30 ) 



Lemma 6.1. The function Gn,p is harmonic on Bi(0). 

Gn,p is harmonic on -Bi(O) away from the coordinate planes. 



4.1 



Proof. By Theorem 
Therefore, it is enough to show that Gn,p is analytic on Bi(0). 



Consider first p = (0,0,0). Then (64) holds on -Bi(O). Since s\ ^ a on £>i(0), 



Gn,p is analytic on Bi(0) \ (A\ U A%) as a composition of analytic functions. In order 
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to show that Gn,p is also analytic at the points of A x U A 2 , one may refer to a classical 
result on "singular curves" of harmonic functions; see Kellogg (1967) (6j Theorem XIII, 
page 271] but we will argue more directly. Since A\ and A 2 are disjoint sets, it is 
clear that B^] (^j is analytic at every point in -Bi(O) fl A\. In order to show that 
^i*n,p( s i)-^2*n,p( S2 ) is analytic at (x',y',z') G -Bi(O) fl A\, we argue as follows. We 
may assume that there is an analytic function w : (ai,a 3 ) — > R such that fif] 
and ^np are restrictions of this function to (ai,a 2 ) and (02,03), respectively. Now 
(si — ai) + (s 2 — Oi) and (si — ai )(s 2 — ai) are analytic functions of (x,y,z) in a 
neighborhood of (x',y',z'). Lemma p^2 implies that E[^ p {si)E^ p {s 2 ) as a function 
of (x,y,z) is analytic at (x',y',z'). It follows that Gn,p is analytic at every point in 
Bi(0) DAi. In the same way, we show that G^j, is analytic at every point in -Bi(O) C\A 2 . 
If p = (0, 0, 1) then we introduce the function 

V a 3 - s 3 if z > 

v : = 

— V a 3 ~ s 3 otherwise. 



It follows from (19), (20) that x is analytic on R 3 \ A%. Then 

G^(x, y, z) = (x 2 + y 2 + z 2 + l)' 1 / 2 'E[]lj Sl )Ei]l p (s 2 )x(x, y, z)w 3 (s 3 ) 

on .Bi(O), where w 3 is analytic at S3 = 03. We then argue as above. 

The other parity vectors p are treated similarly. □ 

Lemma 6.2. Let f : (B e ) 2 — > C, B e = {s G C : \s\ < e}, be an analytic function which 
is symmetric: f(s,t) = f(t, s). Let g, h : (Bs) 3 — >■ B € be functions such that g + h and 
gh are analytic. Then the function f(g(x,y,z),h(x,y,z)) is analytic on (Bs) 3 . 



Substituting tj = fi(sj), j = 2,3, the Hilbert space Hi from Section |5j transforms 
to the Hilbert space Hi consisting of measurable functions g : (ai,a 2 ) x (025^3) - 1 C 
for which 

raz ra,2 s ^ _ 

\\g\\ 2 ■= / / \g(s 2 ,s 3 )\ 2 ds 2 ds 3 < 00. (66) 



ci2 J ai 



u(s 2 )u(s 3 ) 



By Theorem 5.3, for g G Hi and fixed p, we have the Fourier expansion 



g(s 2 ,s 3 ) ^c n)P < p ( S2 )£ji p ( S3 ), (67) 

n 

where the Fourier coefficients are given by 

c n, P =/ f f% f , g(s 2 ,s 3 )E { 2 ] ] (s 2 )E^ ] (s 3 )ds 2 ds 3 . (68) 

Ja 2 Jai UJ {S 2 )U{S 3 ) 

Theorem 6.3. Consider the region Di defined by MJfy for some fixed di G (ao, 01). Let 
e be a function defined on its boundary dDi of parity p G {0, l} 3 , and let g(s 2 ,s 3 ) be 
the representation of 

f{x, y, z) := (x 2 + y 2 + z 2 + l) 1/2 e(x, y, z) (69) 
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in 5-cyclide coordinates for (x, y, z) G dDi with x,y,z > 0. Suppose g G Hi and expand 



g in the series (61). Then the function u(x,y,z) given by 

n -^l,n,p(«l) 



(70) 



harmonic in Di and assumes the values e on the boundary of Di in the L 2 -sense 
explained below. 



Proof. Let di < d < ai and Si G [d, ai]. Using Theorems 5.4, 5.6 



we estimate 



(i) 



E i%, P (di] 



< \c n , p \Cr n *+ n *B(nl + nl + l), 



where the constants B, C > and r G (0,1) are independent of n and Si G [d, ai], 
s 2 G [a 1; a 2 ], s 3 G [a 2 ,a 3 ]. Since c niP is a bounded double sequence, this proves that 
the series in (70) is absolutely and uniformly convergent on compact subsets of D\. 
Consequently, by Lemma 6.1| u(x,y,z) is harmonic in Di. If we consider u for fixed 
si G (di,ai) and compute the norm \\u — e\\ in the Hilbert space Hi by the Parseval 
equality, we obtain 



\u 



< 



i cn .pi 



It is easy to see that the right-hand side converges to as si — > d\. Taking into account 
that e and u have the same parity, it follows that u assumes the boundary values e in 
this L 2 -sense. □ 



If e is a function on dDi without parity, we write the function / from (69) as a 
sum of eight functions 

p 

where f p is of parity p. Then the solution of the corresponding Dirichlet problem is 
given by 



where 



u(x,y,z) = {l) a ' P s G<U p ( x ,y,z), 

,p E l,n,p{dl 
S3 - S 2 



n,p 



a 3 ra 2 



g p {s 2 , s 3 )E 2 : ]l yP {s2)E ( 3 1 J i p (s 3 ) ds 2 ds 



(1) 



(71) 



(72) 



i ai u(s 2 )u(s 3 y 

and g p (s 2 ,s 3 ) is the representation of f p in 5-cyclide coordinates. 

We may write the coefficient integral over the surface dDi itself. The 

surface element is dS = h 2 h 3 ds 2 ds 3 with the scale factors h 2 ,h 3 given in (24), (25). 
Using 



hoh 



2'*3 



1 



hi 



x 2 + y 2 + z 2 + 1) 



w Si 



u(s 2 u(s 3 



yS 3 - S 2 
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we obtain from ( 72 ) 



2u(d 1 )E[ 1 l(d 1 ) Jan, hi 



G^l dS, 



n,p 



(73) 



where 



h\ 



16 



[x 2 + y 2 + z 2 - I) 
(d 1 - a ) 2 



+ 



4x 2 



+ 



Ay 2 



Az 2 



(d x - ai) 2 ((ii - a 2 ) 2 (d x - a 3 ) 2 



7. Second two-parameter Sturm-Liouville problem 

We treat the two-parameter eigenvalue problem that appears when we wish to solve the 
Dirichlet problem in ring cyclides. It is quite similar to the one considered in Section 
[5j however, there are also some interesting differences. Consider equation (26) on the 
intervals (ao,ai) and (02,03). We obtain two Sturm-Liouville equations involving two 
parameters 

1 f 3 



(uj(s 3 )w' 3 y - 



ui(si) \16 



u(s 3 ) \1Q 



s 1 + A1S1 + A 2 1 Wx = 0, a < si < at, 



si + A1S3 + A 2 )w 3 = 0, a 2 < s 3 < a 3 . 



(74) 



(75) 



We again simplify by substituting tj = Q(-Sj), Uj(tj) = Wj(sj). Then (74), (75) become 

3 
16 
3 



:{0(tl)r + A 1 0(ti)+A 2 Ui = O 



u 3 - 



16 



{0(t 3 )r + Ai0(t 3 ) + A 2 U 3 = 



bo < h < 61, 
b 2 < t 3 < b 3 . 



(76) 
(77) 



We add boundary conditions 

u'tibo) = ui(6i) = u' 3 (b 2 ) = u' 3 (b 3 ) 



0. 



(7f 



Differential equations (76), (77) together with boundary conditions (78) pose a two- 



parameter Sturm-Liouville eigenvalue problem. In contrast to Section |5j we now have a 
uniformly right-definite problem: 



<f)(t 3 ) - 0(*a) > a 2 - a a > for b < t x < h < t 3 < b 3 



<f>(h) 1 
<f>(h) 1 

We again have Klein's oscillation theorem. 

Theorem 7.1. For every n = (ni,n 3 ) G Nq ; there exists a uniquely determined 
eigenvalue (Ai jn , A 2 _ n ) G R 2 admitting an eigenfunction u\ with exactly n\ zeros in 
(bo,bi) and an eigenfunction u 3 with exactly n 3 zeros in (b 2 ,b 3 ). 

We state a result on the distribution of eigenvalues. 
Theorem 7.2. There are constants A±, A 2 , A 3 > such that, for all n G Nq ; 

-Atinl + 1) < Ai,„ < A 2 (nl + l), (79) 
|A 2j „| < A 3 (n\ + n 2 3 + l). (80) 
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Proof. We abbreviate Xj = Xj jXl . Arguing as in the proof of Theorem 5.2 there are 
t\ G [b , bi] and t 3 G [62, 63] such that 



16 

3 

16 



M*i)} 2 + Ai0(*a) + A 2 
{^ 3 )} 2 + Ai0(* 3 ) + A 2 



2^2 



7T n 



(6l 



2^2 



7T W 



(63 - h 2, 



2 ' 



(81) 
(82) 



By subtracting (81) from (82), we obtain 
3 



2^2 



16 



(W3)} 2 - W1)} 2 ) + Ai(0(t 3 ) - 0(*O) 



7r n 



2^2 



7r n 



(61 - 60) 



M 2 



which implies (79). Now (80) follows from (79) and (81). 



□ 



Let Mi in and u 3>n denote eigenfunctions corresponding to the eigenvalue (Ai jn , A 2 , n )- 
The system of products Ui tTl (ti)u 3 ^(t 3 ), n G N 2 ,, is orthogonal in the Hilbert space H 2 
consisting of measurable functions / : (bo, bi) x (6 2 , 63) — > C satisfying 

(0(t 3 ) - |/(ti,t 3 )| 2 dtidta < 00 



'62 J bo 

with inner product 

fb 3 pbi 



(0(t 3 ) - t 3 )g(ti,t 3 ) dh dt 3 . 



'62 </ bo 

We normalize the eigenfunctions so that 

f&3 ff>l 



(0(t 3 ) - 0(tl)) {«l,n(*l)} Kn(* 3 )} dtl d*3 = I- 

'62 -'feo 

We have the following completeness theorem. 
Theorem 7.3. The double sequence of functions 

Wl,n(*l)«3,n(*3), n € Nq, 

forms an orthonormal basis in the Hilbert space H 2 . 



(83) 



The normalization (83) leads to a bound on the values of eigenfunctions. Since we 



have uniform right-definiteness, the proof is simpler than the proof of Theorem |5.4| 

Theorem 7.4. There is a constant B > such that, for all n G Nq and all t\ G [bo, bi], 
t 3 G [62,63], 

Kn(£iKn(t 3 )| < B(n\ + n\ + 1). 
Let «2, n be the solution of 
3 



16 



{(t>{t 2 )Y + Ai,„0(t 2 ) + A 2 ,n 1*2 = 0, 61 < £ 2 < 62 



J4) 



determined by initial conditions 
tt 2 (&i) = 1, tt' 2 (6i) 
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Theorem 7.5. We have W2,n(*2) > for all t 2 G [bi,b 2 ]. If b\ < C\ < c 2 < b 2 , then 
there are constants C > and < r < 1 such that, for all n G Nq and t 2 G [b\, c-\\, 

u 2 ,n(t 2 ) 



«2,n(C 2 



< Cr 



ni+n 3 



Proof. We abbreviate w 2 = w 2 , n and Aj = Aj jn . We write (84) in the form 
u 2 + Q((f>(t 2 ))u 2 = 0, t 2 e[h,b 2 ], 

where Q is given by (50). According to (81), (82), there are si G (a , a\) and S3 G (a 2 , 03] 
such that 



Q(si) 



n 2 n\ 



2^2 



7r n 



b^ 2 ' 



(&i - 6 ) 2 ' k (63 - "2J 

If s G [si,s 3 ], then Q(s) < L(s), where L(s) is the linear function with L(sj) = Q(sj), 
j — 1, 3. It follows that 

W(* 2 )) < -C(ni + n 3 ) 2 for t 2 G [c a , c 2 ]. (85) 

Again, we apply Lemma 5.7 to complete the proof. □ 

The results of this section remain valid for other boundary conditions. This time 
we will need sixteen sets of boundary conditions labeled by p = (po,Pi,p 2 ,P3) G {0, l} 4 . 



(86) 



These boundary conditions are 














«i(6o) 


= 


if po = 


= 0, 


Ul(6 Q ) = 


= 


if 


Po = 


= 1 


u[(h) 


= 


if pi = 


= 0, 




= 


if 


Pi = 


= 1 


u' 3 (b 2 ) 


= 


if p 2 = 


= 0, 


u 3 (b 2 ) -- 


= 


if 


P2 = 


= 1 


u' 3 (b 3 ) 


= 


if P3 = 


= 0, 


u 3 (h) -- 


= 


if 


P3 = 


= 1 


The initial conditions 


for 


u 2 are 















m 2 (6i) = 1, w' 2 (6i) = if Pl = 0, u 2 {b x ) = 0, u'M = 1 if Pl = 1. (87) 



We denote the corresponding eigenvalues by (A^ „ , A 



denoted by E 



(2) f 



,n,p5 /v 2,n,p; 



The eigenfunctions will be 



1,2,3. 

,(2) 



Summarizing, for i = 1,2, 3, ^ n J p is a solution of (26) on (aj_i, a*) with (Ai, A 2 ) = 
(^i 2 n, P 5 ^2 2 n, P )- The solution (sj has exponent |p at cto; exponent |pi at a 1; and 
it has ni zeros in (a ,ai). The solution E^'y^) has exponent ~pi at a\, and its has 
no zeros in (ai,a 2 ). The solution E^ p (ss) has exponent ^p 2 at a 2 , exponent \p 3 at 03, 
and it has n 3 zeros in (a 2 ,a 3 ). 
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Figure 5. Coordinate surfaces S2.3 = const for a, = i. 



8. Second Dirichlet problem 

Consider the coordinate surface (21) for fixed s = d 2 G (01,02). See Figure [5|a,b,c) 
for a graphical depiction of the shape of this surface. If (x f , y', z') G §2 then the ray 
(x,y,z) = t(x',y',z') , t > 0, is tangent to the surface if and only if (x',y',z') is on the 
surface. If (x',y',z r ) is in the elliptical cone 
Ax 12 Ay' 2 Az' 2 

+ 1 — + - > 0, 



d 2 — 0\ d 2 — CL 2 d 2 — 0-3 

then the ray does not intersect the surface. Otherwise we have two intersections t = ti, t 2 
and t\t 2 = 1. It follows from these considerations that S2 = d 2 describes a connected 
surface of genus 1. The region interior to this surface is 

D 2 = {(x,y,z)eR 3 :s 2 <d 2 }, (88) 

or, equivalently, 

$2 — ^2 — CX^ Qj 2 — Qj 2 (1 2 — O3 

In this section we solve the Dirichlet problem for harmonic functions in D 2 by the 
method of separation of variables. 

i,n,p 



Let p = (po,Pi,P2,P3) G {0, l} 4 and n = (711,713) G N§. Using the functions E, 



introduced in Section [7] we define the internal 5-cyclidic harmonic of the second kind 
G£> p (x,y,z) = (x 2 + y 2 + z 2 + l)- 1 ^E^{a 1 )E^{a 2 )ESi p (s i ) (89) 
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(x 2 + y 2 + , 2 + l) 1 /^) 



to R 3 as a function of parity p. We call a function / of parity p = (po,Pi,P2,P3) if 



f(<Ti(x,y,z)) = (-l) Pi f{x,y,z) 



for i = 0, 1, 2, 3, 



(90) 



using inversion (29) and reflections (30). 



We omit the proof of the following lemma which is similar to the proof of Lemma 



Lemma 8.1. The function Gn}p is harmonic at all points (x, y, z) G R 3 at which s 2 7^ Ci2>' 
see (FM). 



(2) 

Note that S2 < d 2 < «2 in D 2 . Therefore, Gn,p is harmonic in an open set containing 
the closure of D 2 . Geometrically speaking, the set s 2 = a 2 consists of the part of the 
plane y = "outside" the two closed curves in Figure |3j The ring cyclide D 2 passes 
through the y = plane inside those two closed curves. 

Substituting tj = £l(sj), j = 1,3, the Hilbert space H 2 from Section [7] transforms 
to the Hilbert space H 2 consisting of measurable functions g : (ao, a\) x (02,03) — > C 
for which 

CO3 rai 



\9\ 



S3 - si 



\g(s 1 ,s 3 )\ dsids 3 < 00. 



By Theorem 



7.3 



a 2 Ja W(S1)W(S 3 ) 

for g E H 2 and fixed p, we have the Fourier expansion 



(91) 



(92) 



where the Fourier coefficients are given by 



a 3 pax 



,(2) 



^(^L) ^ 1 ' S 3)^l 2 n, p (si) J B3^ p (s 3 ) dSi rfs 3 . 



Theorem 8.2. Consider the region D 2 defined by (88) for some fixed d 2 G (ax, 02)- ^ei 
e be a function defined on its boundary dD 2 , and set 

f{x, y, z) := (x 2 + y 2 + z 2 + l) 1 / 2 e(a;, y, z). (93) 

Suppose that f has parity p G {0, 1} 4 ; and zis representation g(s\,s 3 ) in 5-cyclide 



coordinates is in H 2 . Expand g in the series (92). Then the function u(x,y,z) given by 

-1(2) 7TT^n,p( a ') 2/) z ) 



u(x,y,z) = 



(94) 



is harmonic in D 2 and assumes the values e on the boundary of D 2 in the L 2 -sense 
explained below. 



Proof. The proof is similar to the proof of Theorem 6.3 It uses Theorems 7.4 and 7.5 to 
show that the series in (94) is absolutely and uniformly convergent on compact subsets 
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of D 2 . Consequently, by Lemma 8.1, u(x,y, z) is harmonic in D 2 . If we consider u for 



fixed S2 € (ai, d 2 ) and compute the norm \\u — e|| in the Hilbert space if 2 , we obtain 



ll" ^ V" 1 1° / 1 (fi 2 ) 

|« - e|| < 2^ l c n, P l 1 - 



The right-hand side converges to as s 2 — > d 2 . Hence u assumes the boundary values e 
in this L 2 -sense. □ 

If / is a function without parity, we write / as a sum of sixteen functions 

/= E /p. 

pe{o,i} 4 

where / p is of parity p. Then the solution of the corresponding Dirichlet problem is 
given by 



u{x, y,z) = Y^ ( ^ G^ p {x, y, z) 

n,p -^2,n,p(«2) 



(95) 



where 



n >p 



13 rai 



S3 - Si 



<7 P ( S1 , s 3 )E[ 2 l p ( Sl )E^ p (s 3 ) d Sl ds 3 (96) 



,(2) 



and gp(si,s 3 ) is the representation of / p in 5-cyclide coordinates. We may also write 
c n , P as a surface integral 

1 f e 



where 



nP Md2)Eil P (d 2 )J 9D2 h 2 G »* dS > 
1 f(x 2 + y 2 + z 2 - l) 2 



(97) 



16 



4x 2 4u 
+ ^ + 



2 4z 2 
+ 



(d 2 - a ) 2 (d 2 - ai) 2 (d 2 - a 2 ) 2 (d 2 - a 3 ) 2 



9. Third two-parameter Sturm-Liouville problem 



If we write (26) on the intervals (ao, Oi) and (01,02) in formally self-adjoint form, we 
obtain 

{y){s\)vJ^)' ]— [ ^-s\ + A1S1 + A 2 J wi = 0, a < si < ai, (98) 



u(si) V 16 



(w(s 2 )u? 2 )' + 



s 2 + Xis 2 + \ 2 ) w 2 = 0, ai < s 2 < a 2 . (99) 



u)(s 2 ) \16 

We simplify the equations by substituting t,- = 0(s 3 -), = Wj(sj), where fi(s) is the 



elliptic integral (|42|). Then (|98j), (99) become 

3 



u. 
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{0(ti)} 2 + Ai0(ti) + A 2 Ui = 0, 6 < *i < 61 



«2 + ( 7^ W2)} 2 + Ai0(t 2 ) + A 2 ) u 2 = 0, 6i < t 2 < 62. 



(100) 
(101) 
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Of course, this system is very similar to the one considered in Section [5j Therefore, 
we will be brief in this section. For a given p = (po,Pi,p 2 ) G {0, l} 3 we consider the 
boundary conditions 

u'i(bo) = if p = 0, ui(b ) =0 if p = 1, 

ui(6 a ) = «' 2 (6i) = if Pl = 0, Ul (6i) = U2(6i) = if Pl = 1, (102) 
m' 2 (6 2 ) = if p 2 = 0, u 2 {b 2 ) =0 if p 2 = 1. 

The initial conditions for u 3 are 

u 3 {b 2 ) = 1, u' 3 {b 2 ) = if p 2 = 0, n 3 (6 2 ) = 0, u' 3 (b 2 ) = 1 if p2 = 1. (103) 

We denote the corresponding eigenvalues by (A[ 3 „ p , A 2 3 „ p ), where n = (nx,n 2 ) G Nq. 
The eigenfunctions will be denoted by p (sj) = «j, n (*i), i = 1, 2, 3. 



Summarizing, for i — 1,2, 3, is a solution of (|26|) on (aj_i, a;) with (Ai, A 2 ) 



(Ai 3 n, P 5 A 2 3 n P ). The solution £'ifnp( s i) nas exponent |p at a o, exponent \p\ at ai, and 
it has ni zeros in (a ,ai). The solution -E , 2 3 n >p (s 2 ) has exponent \p\ at ai, exponent \p 2 
at a 2 , and it has n 2 zeros in (ai,a 2 ). The solution £?f] (s 3 ) has exponent |p 2 at a 2 , 
and it has no zeros in (a 2 , 03). 

10. Third Dirichlet problem 

Consider the coordinate surface (21) for fixed s = d 3 G (a 2 ,a 3 ). See Figure [5|d,e,f) 
for a graphical depiction of the shape of this surface. If (x' , y', z') G §2 then the ray 
(x, y, z) = t(x', y' , z') , t > 0, is tangent to the surface if and only if (x', y 1 ', z') is on the 
surface. If (x',y',z') is in the elliptical cone 

Ax' 2 Ay' 2 Az' 2 

+ ~ — + < 0, 



d 2 — ax d 2 — a 2 d 2 — a 3 
then the ray intersects the surface twice at t — t\,t 2 with t\t 2 = 1. Otherwise, there is 
no intersection. Therefore, the coordinate surface s 3 = d 3 consists of two disjoint closed 
surfaces of genus separated by the plane z — 0, and they are mirror images of each 
other under the reflection a 3 . 

We consider the region inside the surface S3 = d 3 with z > 

D 3 = {(x,y,z) : z> 0, s 3 < d 3 }, (104) 

or, equivalently, 

C3 = {(*.^):z>0 ^ 2 + f + * 2 - 1) V -^ + - 7 ^ + -^. 

Ct 3 — CLq Ct 3 — CL\ CL 3 — CL 2 Ct 3 — CL 3 

Next, we solve the Dirichlet problem for harmonic functions in D 3 by the method of 
separation of variables. 

2 TTciT-,rr fV,^ fnnnfmno Z?( 3 ) 



Let p = (po,Pi,p 2 ) G {0, l} 3 and n = {rii,n 2 ) G Nq. Using the functions E, : 



i,n,p 



introduced in Section [9] we define the internal 5-cyclidic harmonic of the third kind 

G®(x,y,z) = (x 2 + y 2 + z 2 + ir^E^is^E^Js^ljs,) (105) 
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for (x,y,z) G -Bi(O) with x,y,z> 0. We extend the function 

(x 2 + y 2 + z 2 + l^ 2 G^ p (x, y ,z) 

to the half-space {(x, y, z) : z > 0} as a function of parity p. We call a function / of 
parity p = (p ,Pi,P2) if 

/Mar.y,*)) = (-l)«/(ar,y,;z), for z = 0,1, 2 (106) 

using the inversion <To and the reflections a%, a 2 . As before we have the following lemma. 

Lemma 10.1. The function Gn,p is harmonic on {{x, y,z) : z > 0). 

We have the Hilbert space H 3 consisting of measurable functions g : (a ,cti) x 
(d, a 2 ) — > C for which 

ra-2 ra.i _ 

\\g\\ 2 := / / —7 ^— — -^r \g(s 1 ,s 2 )\ 2 ds!ds 2 < 00. (107) 
For g G H3 and fixed p, we have the Fourier expansion 

g(s 1 , 82 ) ~^c^ p e£I p ( 8i )4 3 I p (s 2 ), (108) 

n 

where the Fourier coefficients are given by 

c D)P = / / 2 1 g(s x , s 2 )E[*l (si)E%l (s 2 ) ds l ds 2 . (109) 

Theorem 10.2. Consider the region D 3 defined by (104) for some fixed d% G (02,03). 
Let e be a function defined on its boundary dD 3 , and set 

f(x, y, z) := (x 2 + y 2 + z 2 + l) 1/2 e(x, y, z). (110) 

Suppose that f has parity p G {0, l} 3 , and its representation g(si,s 2 ) in 5-cyclide 
coordinates is in H 3 . Expand g in the series (108). Then the function u(x,y,z) given 
by 

u(x, y,z) = J2 Jr tJ - G%{x, y, z) (111) 

is harmonic in D3 and assumes the values e on the boundary of D3 in an L 2 -sense. 
If / is a function without parity, we write / as a sum of eight functions 

P e{o,i} 3 

where f p is of parity p. Then the solution of the corresponding Dirichlet problem is 
given by 

u(x, y,z) = J2 -^— G^(x, y, z), (112) 

n,p ^3,n,p("3j 

where 

c n, P = / / -^\-^9p(suS 2 )E ( ^l p (s 1 )E^l p (s 2 )ds 1 ds2 (113) 

Jai Ja W{Si)U){S2) 
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and g p (si, S2) is the representation of f p in 5-cyclide coordinates. Alternatively, we have 
c n ,p = ^ / ^G^dS, (114) 

where 

2 _ 1 / {x 2 + y 2 + z 2 -l) 2 4x 2 4y 2 4z 2 \ 

3 " 16 V (d 3 - a ) 2 + (d 3 - Gl ) 2 + (d 3 - a 2 ) 2 + (rf3-a 3 ) 2 y' ' 
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